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We investigate the tunable vibration filtering properties of one-dimensional diatomic 
granular crystals composed of arrays of stainless steel spheres and cylinders interact- 
ing via Hertzian contact. The arrays consist of periodically repeated three-particle 
unit cells (steel-cylinder-sphere) in which the length of the cylinder is varied system- 
atically. We apply static compression to linearize the dynamic response of the crystals 
and characterize their linear frequency spectrum. We find good agreement between 
theoretical dispersion relation analysis (for infinite systems), state-space analysis (for 
finite systems), and experiments. We report the observation of up to three distinct 
pass bands and two finite band gaps and show their tunability for variations in cylin- 
der length and static compression. 

PACS numbers: *43.58.Kr, 45.70.-n, 46.40. Cd 
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I. INTRODUCTION 



The presence of band gaps, a characteristic of wave propagation in periodic structures, 
has been studied in a wide array of settings involving phononic crystals, photonics, and 
plasmonicsi - -. Materials exhibiting band gaps are of particular interest as they forbid and 
allow the propagation of waves in selected frequency ranges (pass and stop bands), and in the 
case of elastic wave propagation (in composites or multilayered structures) have previously 
been proposed for use in acoustic filters, vibration isolation applications, and rectification 
of acoustic energy flux- - -. 

Chains composed of elastic particles in close contact with each other, or "granular crys- 
tals", have gained much attention with respect to elastic wave propagation in nonlinear 
media. The nonlinearity in granular crystals results from the Hertzian contact between two 
elastic spherical (or spherical and cylindrical) particles in compression and from a zero tensile 
strength. The contact stiffness is defined by the geometry and material properties of the par- 
ticles in contact^. In this type of system the level of nonlinearity present can be controlled 
by the amount of static compression applied to the chain, resulting in a dynamic response 
which encompasses the linear, weakly nonlinear, and strongly nonlinear regimes". This 
simple means of controlling their dynamic response has made granular cystals a perfect test 
bed for the study of nonlinear phenomena, including the emergence of coherent structures 
such as solitary waves", discrete breathers", shock waves', and linear /nonlinear defect 
modes". Additionally, granular crystals have been shown to be useful in application to en- 
gineering endeavors including shock and energy absorbing layers""— , actuating devices', 
acoustic lenses' and sound scramblers". 

Previous studies involving statically compressed granular crystals composed of one- 
dimensional (ID) periodic (monoatomic and diatomic with a two particle unit cell) arrays 
of glueo~, welded', and elastically compressed spherical particles'" have been shown 
to exhibit tunable frequency vibrational band gaps. In this manuscript we study statically 
compressed ID diatomic granular crystals composed of periodic arrays of stainless steel 
sphere- cylinder- sphere unit cells. We employ theoretical models to estimate the dispersion 
relation of the crystals, we numerically validate their dynamic response using state-space 
analysis and verify experimentally the crystal's acoustic transmission spectrum. For such 
configurations we experimentally report the presence of a third distinct pass band and a 
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second finite band gap. We show tunability and customization of the response for variation 
of the cylinder length and static compression. 

II. EXPERIMENTAL SETUP 

We assemble five different ID diatomic granular crystals composed of three-particle, 
sphere-cylinder-sphere, repeating unit cells as shown in Fig. Ufa). The chains are 21 par- 
ticles (7 unit cells) long. The particles (spheres and cylinders) are made from 440C stain- 
less steel, with radius R = 9.53 mm, elastic modulus E = 200 GPa, and Poisson's ra- 
tio v = 0.3 — . Each of the five chains is assembled with cylinders of a different length, 
L = [9.4, 12.5, 15.8, 18.7, 21.9] mm. The mass of the spherical particles is measured 
to be m = 27.8 g and the mass of the cylindrical particles is measured to be M = 
[20.5, 27.3, 34.1, 40.7, 47.8] g for each of the corresponding cylinder lengths. 

We align the spheres and cylinders in a horizontal ID configuration using a containment 
structure of four polycarbonate rods (12.7 mm diameter). We hold the polycarbonate rods 
in place with polycarbonate guide plates spaced at intervals of 1 unit cell. We apply low 
amplitude broadband noise to the granular crystals using a piezoelectric actuator mounted 
on a steel cube of height 88.9 mm which is fixed to the table. We visualize the evolution of 
the force-time history of the propagating excitations using a calibrated dynamic force sensor. 
The force sensor is composed of a piezo-electric disk embedded with epoxy inside two halves 
of a R = 9.53 mm, 316 stainless steel sphere (of elastic modulus 193 GPa, and a Poisson 
ratio of 0.3—). The sensor is constructed so as to approximate the mass, shape, and contact 
properties of the spherical particles in the rest of the crystat^i^^. rp^ e g^gg^^g^ f orce 
sensor is calibrated against a commercial dynamic force sensor, and has a measured total 
mass and resonant frequency of 28.0 g and 80 kHz, respectively. We insert the dynamic 
force sensor in place of the last particle, located at the opposite end of the crystal from the 
actuator. We condition its output with a 30 kHz cutoff 8-pole butterworth low pass filter 
and voltage amplifier. 

At the opposite end of the crystal with respect to the piezoelectric actuator, we apply a 
static compressive force, F , using a soft (compared to the contact stiffness of the particles) 
stainless steel linear compression spring (stiffness 1.24 kN/m). In this case we can approx- 
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imate this boundary as a free boundary. The static compressive force applied to the chain 
is adjusted by positioning, and fixing to the table, a movable steel cube of height 76.2 mm 
so that the soft linear spring is compresssed. The resulting applied static load is measured 
with a static load cell placed in between the steel cube and the spring. 



FIG. 1. [Color online] a) Schematic of experimental setup, b) Schematic of the linearized model 
of the experimental setup. 

III. THEORETICAL SETUP 

A. Dispersion Relation 

We model a ID diatomic crystal composed of n sphere-cylinder-sphere unit cells (and N 
particles) chain of nonlinear oscillators^: 

mini = ai-ij [6i-i } i + ui-x - ui] v + 



where [Y} + denotes the positive part of Y; the bracket takes the value Y if Y > 0, and 
if Y < 0. This represents the tensionless characteristic of our system: when adjacent 
particles are not in contact, there is no force between them. The above model assumes that 
the particles act as point masses. This is valid as long as the frequencies of the applied 
vibrations are much lower than the frequencies of the natural vibrational modes of the 
individual particles^. Here, ui is the displacement of the Zth particle around the static 
equilibrium, Si-ij is the static overlap between the (/ — l)th and the Zth particles, and mi 
is the mass of the Zth particle (where I is the index of the Zth particle in the chain counted 



a) 




Linearized sphere-sphere & cylinder-sphere contacts 



(1) 



<xi,i+i[8i,i+i +u t - u i+1 ] p + , 



4 



from the piezoelectric actuator end, and I 6 {1, • • • , 3n}). As per Hertz's contact law, the 
coefficients a depend on the geometry and material properties of the adjacent particles and 
on the exponent p (here p = 3/2) ~. Here, in the case of the sphere- cylinder- sphere unit 
cell, we need to account for two different values of the contact coefficients a, corresponding 
to the sphere- cylinder and the sphere-sphere contacts, where: 

- A - (?) 

& sphere, cylinder cylinder, sphere „/-. 2\ ' \ / 



4 - EV2R , v 

(-^sphere, sphere -^2 o\ " \ / 

3(1 — Z^J 

For this case, it can be seen that A\ = V2A2. Furthermore, for Hertzian contacts, 
under a static load F Q , we can define the static overlap for the sphere- cylinder con- 
tact aS 'sphere, cylinder — ^cylinder, sphere — 

(F /Ai) 2 ^ 3 , and for the sphere-sphere contact as 
$ sphere,s P here = (F0/A2) 2 / 3 — Considering small amplitude dynamic displacements as 
compared to the static overlap, one can linearize the equations of motion (Eq. [1]). For the 
studied sphere- cylinder- sphere unit cell the particles' linearized equations of motion are: 

mil 3j -2 = p2[u 3 j-3 - W3J-2] - Pl[u 3 j-2 - V>3j-l], 
Muzj-l = Pl[u 3 j-2 - «3i-l] - #L[«3i-l - U 3 j], 

(4) 

mu 3j = /3i[«3i-l - u 3j] ~ A^Kj - u 3j+1 ], 



where j is the number of the jth unit cell (j G {1, • • • , n}), m is the mass of a spherical 
particle, M is the mass of a cylindrical particle, (3i = ^A 2 ^ 3 F^ 3 is the linearized stiffness 

n 9/3 1/3 

between a spherical and cylindrical particle, and fa = \A 2 F is the linearized stiffness 
between two spherical particles. The dispersion relation for a diatomic (two particle unit 
cell) granular crystal is known to contain two branches (acoustic and optical)^. Here we use 
a similar procedure to calculate the dispersion relation for a diatomic crystal with a three 
particle unit cell. 

We substitute the following traveling wave solutions where k is the wave number, u is the 
angular frequency, and a = L + AR - 25 spherejCyHnder - 5 sphere)Sphere is the equilibrium length 
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of the sphere- cylinder- sphere unit cell: 

u 3j _ 2 = Ue i{kaj+UJt) 



(5) 

U 3j = W j(kaj+ut)^ 



into Eqs. @. U, V, and W are the wave amplitudes, and are constructed complex so as 
to contain both the amplitude and phase difference for each particle within the unit cell. 
Solving for a nontrivial solution we obtain the following dispersion relation: 

= - 2/3! 2 /3 2 + ft (ft + 2/3 2 )(2m + M)u 2 

-2m(/3 2 M + /3 1 (m + M))w 4 (6) 

+ m 2 Muj 6 + 2/3 1 2 /3 2 cos(afc). 
In Fig. [2] (a), we plot the dispersion relation (Eq. [6]) for the previously described sphere- 
cylinder-sphere unit cell granular crystal with cylinder length L = 12.5 mm (M = 27.3 g) 
subject to an F = 20 N static load. Three bands of solutions (or propagating frequencies) 
can be seen; the lowest in frequency being the acoustic band, followed by lower and upper 
optical bands. Frequencies in between these bands are said to lie in a band gap (or forbidden 
band). Waves at these frequencies are evanescent, decay exponentially, and cannot propagate 
throughout the crystal 1 . 

If we solve the dispersion relation, Eq. (JSJ), for when k — - and k = we obtain the 
following cutoff frequencies: 



/ c 2 i = o 

fc,2 
p2 



2 = /3i + 2/3 2 
/3i(2m + M) 



47r 2 m 
r i(2m , 

/c ' 3 ' An 2 mM 
Pi 



cA 



2 

c,5 



/, 



+ 



47r 2 m' 

0x(2m + M) + 2/3 2 M (7) 
M 2 mM 

yZ-lGP^mM + (2/3im + ftM + 2/3 2 M) 2 

87T 2 mM 
/?i(2m + M) +2/3 2 M 

V^-lGftftrnM + (2ftm + ftM + 2/3 2 M) 2 
87r 2 mM 
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fc,i, fc,2, and / c>3 correspond to k = and / C)4 , / c , 5 , and / C)6 to A; = f . In Fig. [2] (a), we 
label the six cutoff frequencies (Eqs. [7]) for the previously described granular crystal with 
cylinder length L = 12.5 mm (M = 27.3 g) subject to a Fq = 20 N static load. 

From Eqs. [7J it can be seen that the cutoff frequencies are tunable through the variation 
of particle masses m and M, and the linearized stiffnesses (3\ and 02 (thus tunable with 
changes in geometry, and static compression Fq). In Fig. 2(b) we plot the cutoff frequencies 
in Eqs. ([7]) as a function of cylinder length for fixed F =20 N static compression, and in 
Fig. 2(c) as a function of static compression (F = [20, 25, 30, 35, 40] N) for fixed cylinder 
length L = 12.5 mm (M = 27.3 g). The lines represent the cutoff frequency solutions 
(/ Cj 2 is dashed to clarify the nature of the intersection with / Cj3 , and the shaded areas are 
the pass bands). It can be seen that within our frequency range of interest, two of the 
cutoff frequency solutions coincide at specific cylinder lengths. The intersection between / Cj4 
and / Cj 5 can be found to occur at M/m = 1^ and the intersection between f c ^ and / C) 3 at 
M/m = (2 — |^). Notice, however, that aside from these special parameter values where the 
above intersections occur, the spectrum preserves the three pass bands with two associated 
finite bandgaps between them. 



B. State-space approach 

In addition to the dispersion relation previously calculated for an infinite system, we 
study the finite linearized system corresponding to our experimental setup as shown in 
Fig. Hlb). We model the actuator boundary of our system as a fixed 440C steel wall. We 
model the other end of the chain as a free boundary, as the stiffness of the spring used for 
static compression is much less than the characteristic stiffness of the particles in contact. 
The linearized equations of motion for the finite system are the same as Eqs. fl4]), except the 
equations for the first and last particles which are given by the following expressions: 

mu x =F\ — - 0\\u\ - u 2 ], 

mui\ = /?i[«20 - u 2 i\. (8) 
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FIG. 2. a) Dispersion relation for the described sphere-cylinder-sphere granular crystal with cylin- 
der length L = 12.5 mm (M = 27.3 g) subject to a -Fo = 20 N static load. The acoustic branch 
is the dashed line, the lower optical branch is the solid line, and the upper optical branch is the 
dash-dotted line. Cutoff frequencies for granular crystals corresponding to our experimental con- 
figuration (b) varying the length L (and thus mass) of the cylinder with fixed ^0=20 N static 
compression, and (c) varing the static compression (Fq = [20, 25, 30, 35, 40] N) with fixed 
L = 12.5 mm cylinder length (M = 27.3 g). Solid lines represent the six cutoff frequency solutions. 
f C) 2 is dashed to clarify the nature of the intersection with f c ^. Shaded areas are the propagating 
bands. 

where Fi is the force applied to the first particle by the actuator. Next, we apply the 
state-space approach, which can be written as^: 

x = Ax + Bu 

y = Cx + Du, (9) 



S 



where vectors x, u, and y are the state, input, and output vectors, respectively. Matrices A, 
B, C, and D are called state, input, output and direct transmission matrices, respectively. 
We choose as an input to the system the force Fx, i.e., u = Fx and as an output y = F^ = 
ffi[n2o-"2i] ^ a veraged force of the two contacts of the last particle (which is analogous to 
what is measured by the embedded dynamic force sensor in our experimental setup ^L2iMi2£i22. 
Thus, for the linear system of Fig. [T^b), we obtain 
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and D = 0. is a zero matrix and I is the identity matrix (both of size N x N). 

We use the formulation in Eqs. [9] with MATLAB's (R2008b) bode function to compute 
the bode diagram of the frequency response for the experimental configurations described. 
The bode diagram is the magnitude of the transfer function H(s) = D + C(sl — A) _1 B, 



where s 
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. We plot the bode transfer function \H(iuS) \ for the five previously described 
diatomic (three-particle unit cell) chains with varied cylinder length for fixed -F =20 N static 
compression, (Fig. |3]^a)), and with varied static compression (F = [20, 25, 30, 35, 40] N) 
for fixed cylinder length L = 12.5 mm (M = 27.3 g) (Fig. [3]^b)). 

We truncate the visualization in Fig. E] below —40 dB and above 20 dB as a visual 
aid to maintain clarity of the frequency region of interest. This resembles experimental 
conditions, as the noise floor of our measurements is approximately —38 dB (as can be seen 
in Fig. H]) and the presence of dissipation in our experiments reduces the sharpness of the 
resonant peaks in contrast to those predicted by the state-space analysis. Attenuating and 
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propagating frequency regions for this formulation match well with the cutoff frequencies 
of the infinite system (see Eqs. ©), denoted by the solid lines plotted in Fig. [3j The high 
amplitude (bright) peaks correspond to the eigenfrequencies of the system, the modes of 
which are spatially extended. However, for certain cylinder lengths, we also observe an 
eigenfrequency located in the second gap of the linear spectrum (denoted by an arrow in 
Fig. [3(a)). These modes result from the break in periodicity due to the presence of the 
actuator "wall" (acting like a defect in the system). In our setup, it can be seen in Fig. [1(b) 
that the first particle (which is spherical) is coupled to both its nearest neigbors via springs 
characterized by spherical-planar contact (/?i). This is unique within the chain and forms a 
type of locally supported defect mode. When the frequency of this mode lies within a band 
gap the mode becomes spatially localized around the first particle and its amplitude decays 
exponentially into the chain. Furthermore, as our chains are relatively short and the gap 
that the localized modes occupy relatively narrow (in frequency), the spatial profile is found 
to be almost similar to the extended modes. This suggests that it may be experimentally 
difficult to differentiate these modes from their extended counterparts. 




Cylinder Mass (g) Static Load (N) 



FIG. 3. [Color online] Bode transfer function (\H(iuj)\) for the experimental configurations: (a) 
the five diatomic (three-particle unit cell) granular crystals with varied cylinder length for fixed 
Fq=20 N static compression, and (b) the fixed cylinder length L = 12.5 mm (M = 27.3 g) granular 
crystal with varied static load. Solid white lines are the cutoff frequencies calculated from the 
dispersion relation of the infinite system. The black arrow in (a) denotes the eigenfrequency of a 
locally supported defect mode. 
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IV. EXPERIMENTAL LINEAR SPECTRUM 



We experimentally characterize the linear spectrum of the previously described diatomic 
chains with sphere-cylinder-sphere unit cells for varied cylinder length and static load. We 
apply a low-amplitude (approximately 200 mN peak) bandwidth limited (3 — 15 kHz) noise 
excitation with the piezoelectric actuator. We measure the dynamic force using a sensor 
embedded in the last particle of the granular crystal as shown in Fig. [TJ We compute the 
power spectral density (PSD 35 ) of the measured dynamic force history over 1.3 s intervals, 
and average the PSD over 16 acquisitions. We normalize the averaged PSD spectrum by 
the average PSD level in the 3 — 7.5 kHz range of the L = 12.5 mm (M = 27.3 g), F = 
20 N granular crystal response to obtain the transfer functions shown in Fig. [4] and Fig. |5J 
More specifically, Fig. H] shows the experimental transfer function in more detail for the 
L = 12.5 mm (M = 27.3 g), F = 20 N granular crystal. 

As in^, we observe that the experimentally determined spectra are upshifted in frequency 
from the theoretically derived spectra for all configurations tested. Because of this we use 
the measured spectra to extract the effective elastic properties of our system. For the F = 
[20, 25, 30, 35, 40] N, fixed cylinder length L = 12.5 mm (M = 27.3 g) granular crystals, 
we measure the frequencies of the —10 dB level of the PSD transfer function corresponding 
to the second band gap {f c ,2 and f c ,z)- We use these experimentally determined frequencies 
to solve for two experimentally determined Hertzian contact coefficients of our system Ai i£xp 
and A 2t exp using the previously described equations for Ai, A 2 , 0i, 02, fc,2, and f c ^. An 
example of the determination of f Cy2 and / Cj 3 for the L = 12.5 mm (M = 27.3 g), F = 20 N 
granular crystal is shown in Fig. |H We compare the experimentally determined A ltexp 
and A 2 .exp to the theoretically determined A\ and A 2 in Table [J (error ranges indicate the 
standard deviation resulting from the measurements at the five different static loads). As 





A x [N//im 3 / 2 ] 


A 2 [N/mn 3 / 2 ] 


Theory 


14.30 


10.11 


Experiments 


18.04 ±0.44 


11.48 ±0.06 



TABLE I. Hertz contact coefficients derived from standard specifications 32 {A\ and A 2 ) versus 
coefficients derived from the measured frequency cutoffs (^4i,exp and ^42,exp), for the (Fo = 
[20, 25, 30, 35, 40] N) fixed cylinder length L = 12.5 mm (M = 27.3 g) granular crystals. 
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FIG. 4. [Color online] Experimental transfer function for the L = 12.5 mm {M = 27.3 g), -Fo = 20 N 
granular crystal. The horizontal dashed line is the —10 dB level used to experimentally determine 
the f c ,2 and / Cj 3 band edges which are denoted by the vertical dashed lines 

the equations for the five non-zero cutoff frequencies (see Eqs. [7]) in our granular crystals are 
dependent on some combination of A\ and A2, the choice of using f c ^ and f c $ to solve for A\ 
and A2 is not unique and other combinations of cutoff frequencies could be used similarly. 

In previous work 1 ™, numerous possible explanations for the upshift in the spectrum have 
been identified. Uncertainty in the standard values of material parameters^ or deviations 
in the local radius of curvature due to surface roughness could result in the material be- 
having more stiffly 36 . Deviations from Hertzian behavior could potentially be caused by the 
dynamic loading conditions, non-Hookean elastic dynamics (due to nonlinear elasticity or 
plasticity), or dissipative mechanisms (viscoelasticity, solid friction), and could result in a 
shift in the exponent p or in the effective contact coeffcient A— r — . We also observe that the 
contact coefficient A between the cylindrical and spherical particles has the larger deviation 
from theory. This deviation could be attributed mainly to the cylindrical particles, due to 
characteristics not shared by the spherical particles. Such characteristics could include: sur- 
face roughness particular to the manufacturing process of the cylindrical particles, or plastic 
deformation occuring closer to the surface as compared to spherical particles. 

In Fig. [5], we plot the experimentally determined PSD transfer functions for the five previ- 
ously described diatomic (three-particle unit cell) chains with varied cylinder length for fixed 
F =20 N static compression (Fig. Eta)), and static compression F = [20, 25, 30, 35, 40] N, 
for fixed cylinder length L = 12.5 mm (M = 27.3 g) (Fig. E^b)). We plot with solid white 
lines the cutoff frequencies from the dispersion relation calculated using the experimentally 
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determined Hertz contact coefficients A^ exp and A 2) ex P - We observe good agreement between 
the semi-analytically derived cutoffs (i.e, from the theoretical dispersion relation but using 
Ai >exp and A 2exp ) and the experimental spectra. By comparing Fig. E]to Fig. Owe observe 
good qualitiative agreement between the numerical (state-space) and experimental spectra. 
Comparing the experimentally and theoretically determined cutoff frequencies, we observe 
an average (over all experimental configurations) upshift in the experimental frequency cut- 
offs versus the theoretically determined frequency cutoffs of: 5.8% in / Cj2 , 8.1% in / C)3 , 8.1% 
in / C)4 , 5.4% in / c>5 , and 7.0% in f cfi . 

The demonstrated attenuation of the elastic wave propagation in frequency regions lying 
within the band gaps of the granular crystals shows that such systems have potential for 
use in a wide array of vibration filtering applications. Furthermore, the tunability displayed 
(achievable from material selection, shape, size, periodicity and application of static com- 
pression) offers significant potential for attenuating a wide spectrum of undesired frequencies. 





Cylinder Mass (g) Static Load (N) 



FIG. 5. [Color online] Experimental PSD transfer functions for the experimental configurations 
described in Fig. [3l (a) the five diatomic (three-particle unit cell) granular crystals with varied 
cylinder length for fixed Fq=20 N static compression, and (b) the fixed cylinder length L = 12.5 mm 
(M = 27.3 g) granular crystal with varied static load. Solid white lines are the cutoff frequencies 
from the dispersion relation using experimentally determined Hertz contact coefficients Ai i6xp and 

Al^exp- 
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V. CONCLUSIONS 



In this work, we describe the tunable vibration filtering properties of a ID granular 
crystal composed of periodic arrays of three-particle unit cells. The unit cells are assembled 
with elastic beads and cylinders that interact via Hertzian contact. Static compression is 
applied to linearize the dynamics of particles interaction and to tune the frequency ranges 
supported by the crystal. We measure the transfer functions of the crystals using state-space 
analysis and experiments, and we compare the results with the corresponding theoretical 
dispersion relations. Up to three distinct pass bands and three (two finite) band gaps are 
shown to exist for selected particle configurations. The tunability of the band edges in 
the crystal's dispersion relation is demonstrated by varying the applied static load and the 
cylinder length. 

In the present work, we restrict our considerations to the study of near linear, small 
amplitude excitations. A natural extension of this work would involve the examination of 
nonlinear excitations within the bandgaps of such granular chains 14 . In particular, it would 
be relevant to compare the properties of localized nonlinear waveforms in different gaps of 
the linear spectrum. Such studies will be reported in future publications. 
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